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ABSTRACT 

We have analyzed the time variability of the wide-band X-ray spectrum of Vela X-1, the brightest 
wind-fed accreting neutron star, on a short timescale of 2 ks by using Suzaku observations with an 
exposure of 100 ks. During the observation, the object showed strong variability including several 
flares and so-called "low states" , in which the X-ray luminosity decreases by an order of magnitude. 
Although the spectral hardness increases with the X-ray luminosity, the majority of the recorded flares 
do not show any significant changes of circumstcUar absorption. However, a sign of heavy absorption 
was registered immediately before one short flare that showed a significant spectral hardening. In the 
low states, the flux level is modulated with the pulsar spin period, indicating that even at this state the 
accretion flow reaches the close proximity of the neutron star. Phenomenologically, the broad-band 
X-ray spectra, which are integrated over the entire spin phase, are well represented by the "NPEX" 
function (a combination of negative and positive power laws with an exponential cutoff by a common 
folding energy) with a cyclotron resonance scattering feature at 50 keV. Fitting of the data allowed 
us to infer a correlation between the photon index and X-ray luminosity. Finally, the circumstellar 
absorption shows a gradual increase in the orbital phase interval 0.25-0.3, which can be interpreted 
as an impact of a bow shock imposed by the motion of the compact object in the supersonic stellar 
wind. 
Subject headings: accretion, stars: neutron, pulsars: individual (Vela X-1), X-rays: binaries 



1. INTRODUCTION 

Accretion-powered pulsars in high-mass X-ray bina- 
ries (HMXBs) arc the best laboratories for studying the 
accretion in the presence of a strong magnetic field. 
Such binary systems harbor neutron stars with a typ- 
ical magnetic field of ^10^^ G, which accrete matter 
from the surrounding environment, and their gravita- 
tional energy is realized through the bright X-ray ra- 
diation. In spite of intensive studies of the X-ray pul- 
sars since the dawn of X-ray astronomy, their emis- 
sion mechanism still remains unclear due to the extreme 
physical conditions induced by strong gravitation and 
magnetic fields. Given the shape of the X-ray spec- 
trum, which can be described by a hard power law with 
a quasi-exponential cutoff at an energy of 20-30 keV, 
Comptonization by hot electrons in the accretion fiow 
can provide a feasible scen ario for the formation of the 
X-ray spectral continu um (jSunvaev fc TitarchukI 119801 : 
iMakishima et al.]ll999f ). However, the relation between 
the radiation and the physical conditions of the accretion 
flow has not been understood yet. 

In order to confirm the hypothesis of the Comptonized 
emission, it is essential to build more detailed physi- 
cal models of the accreted plasma that radiates X-rays. 
Since the spectral variability refiects changes of the phys- 
ical conditions of the plasma, the variability patterns 
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can contain important information about the accreted 
plasma. We therefore analyzed X-ray archive data of 
HMXB Vela X-1 obtained with the Suzaku X-ray ob- 
servatory to pursue the variability of hard X-rays from 
the central engi ne. The first analysis resu lts of this data 
are presented bv lDoroshenko et al.l (j2011h . Owing to the 
wide-band coverage and the low -background capability 
of Suzaku (jMitsuda et al.l [2007|) . we are able to obtain 
unprecedentedly high-quality spectra with a short inte- 
gration tim e of 2 ks. Furthermore, the Hard X-ray Detec- 
tor (HXD) (jTakahashi et al.ll200n iKokubun et"an[2007[ ) 
onboard Suzaku is the only instrument that is able to 
obtain such a short-period spectrum in the hard X-ray 
band above 10 keV with high signal-to-noise ratios. 

Vela X-1 (4U 0900-40) is the best target for such 
study since it is the brightest wind-fed pulsar, which al- 
ways displays strong time variability of the X-ray flux 
level. In particular, the source displays intensive out- 
bursts, e.g., the flaring behavior in hard X-rays was ob- 
served with INTEGRAL up to an inten sity of ~ 5 Crab at 
20-40 keV (jKrevkenbohm et al.|[2008[) . The binary con- 
sists of a B0.5Ib supergiant HD 77581 of radius 30.0i?Q, 
and a neutron star with a sp in period of 283.5 seconds 
(jMcClintock et al.lll976t iKrev kcnbohm ct al. 2008.) on a 
nearly c ircular orbit with a smal l binary separation of 
53.4i?0 (Ivan Kerkwiik et aI]|1995D an d orbital period of 
8.964 davs (|van Kerkwiik et al.lll995D . Since the B-type 
companion drives a strong stellar wind with a mass-loss 
rate of about 2 x IQ-^MQyr-^ (.Wa tanabc et al 200| ) 
and terminal velocity of 1100 km s~^ (jPrinia et al.ll990l) . 
the neutron star is always embedded deeply in the dense 
stellar wind. This creates perfect conditions for an in- 
tensive accretion, which leads to the formation of the 
bright X-ray emission. Assuming a distance of 1.9 kpc 
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to the object (jSadakane et al.lll985D . the typical X-ray 
luminosity of Vela X-1 is ~ 4 x 10^^ erg s~^. 

Though it is beyond doubt that the accretion rate 
to the neutron star is highly variable in such systems, 
the mechanism that governs accretion is still in de- 
bate. Consequently, the physical origin of the X-ray 
variability has not been studied in detail yet. One 
straightforward, possible explanation for the time vari- 
ability is that inhomogeneous structure of the stellar 
wind leads to variations of the density and velocity of 
matter around the neutron star, resulting in a fluctu- 
ation of the accr etion rate (iKrevkenbohm et aLl 120081 : 
iDucci et all 120091 : iFiirst et al.l l2010(). Another interpre- 
tation is that a magnetospheric barrier regulates the ac - 
cretion (|Grebenev fc SunvaevI [20071 : iBozzo et al.l 120081 ). 
Within this scenario, the origin of the variability was 
discussed in iDoroshenko et al.l ()2011l ). who has analyzed 
"off states" of Vela X-1 based on the Suzaku data (these 
data are also used in this paper). 

The present paper is organized as follows. Section[2]de- 
scribes the Suzaku observation of Vela X-1 and the data 
reduction. The analysis results, which include the time 
variability of the spectrum, are presented in Sections [3] 
and 131 Section [5] discusses the time variability, the struc- 
ture of the wind, and spectral features of the neutron 
star radiation. In Section [6] we summarize our findings 
and conclusions. Implications of the obtained observa- 
tional results to detailed m odeling of the centr al engine 
are presented separately by lOdaka et al.l ()2013D . 

2. OBSERVATIONS AND DATA REDUCTION 

The observation of Vela X-1 with Suzaku was per- 
formed in June 2008 for an exposure time of 100 
ks (Principal Investigator: T. Kallman). Information 
on the observation is summarized in Table [1] Ac- 
cording to the ephemeris obtained with INTEGRAL 
(jKrevkenbohm et al]|2008[ ). the Suzaku observation cov- 
ered the orbital phase interval between 0.175 and 0.310, 
which corresponds to 13.5% of the orbital period of the 
Vela X-1 binary system of 8.964 days. Since the target is 
very bright, all of the three CCDs (XIS 0, 1, 3) were op- 
erated under the 1/4 window option, in which 256x1024 
pixels of the CCD are read out every 2 s. 

We used the Suzaku data of Vela X-1 processed by a 
set of software for the Suzaku processing pipeline ver- 
sion 2.2.8.20, and used the HEASoft software package 
version 6.9 together with a calibration database updated 
on 2010-08-12 for data processing as well as basic anal- 
ysis including spectral fitting. For all the data with 
the XIS and the HXD, we reprocessed unfiltered event 
data by ourselves to apply new calibration information 
on the detector gain and to avoid a problem of the clean- 
sis tool, which usually eliminates flickering events of the 
CCDs but sometimes fails with old versions due to its 
detection algorithm. In the data reprocessing, standard 
event screening criteria were applied for both the XIS 
and the HXD. These standard screening procedures in- 
clude selection by event grades and removal of time in- 
tervals such as the spacecraft's passage of the South At- 
lantic Anomaly (SAA), regions of small geomagnetic cut- 
off rigidity (COR), and those of a low elevation angle 
from the Earth. For the XIS, events for an elevation 
angle larger than 5° above the Earth limb and larger 
than 20° above the sunlit Earth limb are selected. The 



HXD data selection requires an elevation angle larger 
than 5° above the Earth limb and the COR larger than 
6GV. 

In the following sections the observational results ob- 
tained with Suzaku from Vela X-1 and details of the data 
analysis performed are described. Due to a noticeably 
better performance of Suzaku below 10 keV, the analy- 
sis was conducted in two ways: (1) soft- X-ray analysis 
below 10 keV using only the XIS data (g]); and (2) the 
XIS/HXD simultaneous wide-band analysis (ijl]). The 
dedicated soft X-ray analysis, independent of any im- 
pact of the hard X-ray instrument, is also justifled from 
a physical point of view, since the soft-X-ray properties 
of the Vela X-1 system are particularly useful for study- 
ing the relation between the circumstellar medium and 
the activity of the compact object from the viewpoint 
of the origin of the time variability. On the other hand, 
the hard X-ray information is the key to grasp the ra- 
diation mechanism from accretion-powered pulsars. In- 
deed, the spectral features of the emission from accreted 
plasma appear above 10 keV in the hard X-ray band, 
while the soft X-ray continuum has a featureless power- 
law shape. Generally, the X-ray pulsar spectrum has a 
quasi-exponential cutoff at 20-30 keV, and often shows 
absorption-like features due to the cyclotron resonance 
in the hard X-ray band. The spectral analysis was per- 
formed by using a spectral fitting software package Xspec 
version 12.6 included in HEASoft. 



TABLE 1 
Summary of the Suzaku observation of Vela X-1 



Property 


Value 


Obs. ID 


403045010 


Coordinates (RA, Dec) (J2000) 


(135.5377, -40.5514) 


Pointing 


XIS nominal 


Start time 


2008-06-17 04:45:17 


Start time (MJD) 


54634.19811 


Observation duration 


145 ks 


Orbital phase 


0.175 - 0.310 


Exposure time 


104 ks 


XIS clocking mode 


Normal 


XIS burst option 


No 


XIS window option 


1/4 window 


XIS Spaced-row Charge Injection 


On 



3. SOFT X-RAY DATA ANALYSIS 

In this section we present analysis of data obtained 
in the soft X-ray energy band with XIS of the Suzaku 
telescope. The performed analysis includes a study of 
the pileup effect in the data (§ 13. ip . of the temporary 
behavior of the source (? l3.2p . and of spectral properties 



3.1. Preparatory Analysis and Photon Pileup Effects 

Primarily, we have examined the effects of photon 
pileup on the CCDs. The photon pileup effect is im- 
portant in the cases if the observed source is very bright: 
two or more photons can strike the same CCD pixel dur- 
ing a single readout and generate an event regarded as 
a non- X-ray event or a higher energy photon. This re- 
sults in an inaccurate estimate of the spectral properties 
and source flux. The modest angular resolution of the 



i J: it 







100 -j ^1 
h I 
80 

60 

40 

20 

"0 20 40 60 80 100 120 140 

time [ks] 

Fig. 1.— Light curve of Vela X-1 obtained with XIS (0.3-12 
keV). The time bin size is 24 s. 



Suzaku X-ray mirrors and a short exposure of 2 seconds 
in a single readout of the CCD with the 1/4 window op- 
tion do reduce the pileup effect. However, since Vela X-1 
is an extremely bright source, a careful treatment of the 
possible photon pileup is necessary. 

First, we have extracted a light curve from the circu- 
lar region centered on the image peak with a radius of 
120" to check the significance of the pileup effect and 
temporal behavior of the source. The light curve, shown 
in Figure [1] displays strong time variability on top of the 
X-ray pulsation with a period of 284 s. The count rate 
reaches about 100 counts s~^ in periods of the flaring be- 
havior. Since operations of the XIS with the "1/4 win- 
dow" and "no burst" options requires a count rate less 
than 50 counts s~^ to avoid a strong pileup, the photon 
events recorded from the Vela X-1 can be considerably 
affected by the pileup effects. 

A simple standard method to correct the pileup effect 
is to remove the center part from the event extraction 
region for light curves or spectra. In order to find an ap- 
propriate extraction region, we performed preparatory 
spectral analysis of the XIS data. We defined five ex- 
traction regions in a shape of a circle or annulus: (1) 
< r < 120" (no removal of the image center); (2) 
20" <r < 120"; (3) 30" <r < 120"; (4) 45" < r < 120"; 
and (5) 60" < r < 120". Here r denotes the radial dis- 
tance from the image center. From these regions, we 
extracted spectra of the XIS every 2268 s, which corre- 
sponds to eight times the spin period of the X-ray pulsar. 
Background spectra were extracted from source-free re- 
gions. Detector response matrices ("rmf") of the XIS 
were generated by xisrmfgen, and an ancillary response 
function ("arf") was simulated by xissimarfgen for each 
extraction region and each time interval. 

It is well-known that soft X-ray emissions from Vela 
X-1 below 10 keV can be characterized by an absorbed 
power law with a strong iron fluorescence line at 6.4 keV. 
Thus, we modeled the XIS spectra by fitting to 
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+ Fpe(i?)), (1) 



where E, A^h, crabs (£'), A, T are photon energy, equiva- 
lent hydrogen column density of absorbing matter, cross 
section of the photoelectric absorption on hydrogen, nor- 
malization, and power law exponent, respectively. Here, 
we used wahs model in Xspec for the absorption model. 
The iron line was modeled by a Gaussian F-pe{E) with 



a very narrow width at 6.4 keV. In the fitting, A'h, F, 
A and the Gaussian normalization were treated as free 
parameters. It is also known that the emission scattered 
from an ambient accretion fiow or stellar wind sometimes 
results in a strong soft X-ray excess (below ~2 keV) on 
top of the absorbed power-law continuum. To reduce ef- 
fects of the soft excess on spectral fitting, we restricted 
the energy range for fitting to 2.5-10 keV. 

The fitting results for XIS are shown in Figure [2] 
for the whole observation period. The spectra extracted 
from the region that includes the image center (black 
line) show lower fiuxes and a harder spectrum particu- 
larly at epochs of high fiux level, as compared to the 
spectra extracted from other regions. This is a straight- 
forward indication of photon pileup. The effects become 
less significant when the central region excluded from the 
data analysis is increased. We considered that the region 
of 45" < r < 120" (blue lines in the figure) is suitable for 
our spectral and timing analysis given that the pileup ef- 
fects are negligible and this region contains enough events 
to be statistically significant. The analysis of the XIS 
data described in the following sections is based on the 
events extracted from the region of 45" < r < 120" for 
aU the CCDs. 

3.2. Light Curves 

Since the X-ray emission from Vela X-1 displays a 
strong time variability in addition to the X-ray modu- 
lation caused by the pulsar (i.e., that due to the pulsar 
rotation) , it is essential for understanding of the emission 
nature to examine the obtained light curves. In Figure[3l 
we show light curves obtained with the XIS for different 
energy ranges: 0.6-2 keV, 2-5 keV, and 5-10 keV. The 
light curves were extracted by summing up the count 
rates of XIS and XIS 3. The light curve time bin was 
selected to coincide with the spin period of the neutron 
star, i.e. 284 s. Since the background levels are negligi- 
ble and no flaring of the background occurred during the 
observation, no background subtraction is necessary. 

The count rates in the three different energy ranges 
provide us with the spectral evolution of Vela X-1. The 
lowest-energy band (0.6-2 keV) is heavily affected by 
photoelectric absorption in the accretion flow or stellar 
wind, while the higher energy bands are not signiflcantly 
affected by the absorption. In order to characterize the 
spectral properties at each moment, we have introduced 
the hardness ratio between the hard and soft bands as 



hardness ratio = 



h-s _ H/H - S/S 
h + s " H/H + S/S 



(2) 



Here, h and s are count rates normalized to the average 
count rates in the hard and soft bands, respectively; H 
and S are the corresponding rates in these bands, and the 
bar indicates the value averaged over the whole observa- 
tion. This hardness ratio can be used as an indicator 
of the spectral hardness, ranging between —1 and -1-1. 
In addition to the light curves of count rate, we show 
the hardness ratios for two different sets of energies in 
Figure [31 The hardness ratio calculated for the bands 
of 0.6-2 kcV and 2-5 keV is sensitive to the amount of 
absorbing medium, or Nn. On the contrary, the ratio for 
2-5 keV and 5-10 keV is sensitive to the intrinsic spectral 
slope rather than to the absorption impact. 
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Fig. 2. — Best-fit values of the equivalent hydrogen column density Nn and the photon index F obtained from the XIS spectra as 
functions of time for different extracting regions together with calculated flux and reduced x^ of the fit. The parameters from the extraction 
regions of < r < 120" (black), 20" < r < 120" (red), 30" < r < 120" (green), 45" < r < 120" (blue), and 60" <r < 120" (magenta) are 
shown. The error bars with A^h and photon index F show la uncertainties. 
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Fig. 3. — Light curves from XIS for three different energy ranges 
of 0.6-2 keV, 2-5 keV, and 5-10 keV. Hardness ratios between the 
different energy ranges are also presented. The data count rates 
were obtained by summing those of XIS and XIS 3. The time 
bin size is 284 s, the same as the spin period of the neutron star. 



During the 140-ks-long observation of Vela X-1, several 
flares, which were a few times brighter than the average 
flux level, have been detected. Most of them had a dura- 
tion of ^10 ks and the hardness ratio does not show any 
significant spectral changes during these flares. In con- 



trast, the flare at t = 110 ks had different characteristics: 
it has been characterized by a very short (about 1 ks) du- 
ration accompanied by a significant spectral hardening. 
Figure 3] shows light curves around the fiare with short 
time bins of 28 s. It can be seen that the flare around 
t = 110 ks consisted of two shorter sub-flares at t = 109.3 
ks and t — 110.4 ks. The former displayed a significant 
spectral hardening while the hardness ratio of the latter 
was similar to the typical or average value ('^0.0). Inter- 
estingly, a strong attenuation of the soft band, probably 
caused by heavy absorption, was observed just before the 
first sub-flare. 

Besides the flaring behavior, the object sometimes 
dropped into "low states" , in which the count rates were 
about an order of magnitude smaller than the average. 
Three periods of the "low states" are seen in Figure [31 
t =10 ks, 30 ks, and 40 ks. The "low state" observed at 
f = 40 ks lasted for 3.5 ks, while the duration of others 
was about 1 ks. Apparently, a spectral softening was oc- 
curring during all the "low states" . This spectral feature 
might be essential to understand the physical situation 
of X-ray pulsars in "low states" . Figure [5] shows light 
curves around the low state at i = 40 ks with short time 
bins of 28 s. In the low state. Vela X-1 still displays X- 
ray pulsations, clearly indicating that the accretion flow 
reaches near the magne tic poles on the stellar surface. 
IDoroshenko et ahl (j2011|) reported detailed data analysis 
and its interpretation based on the same Suzaku data of 
these low states. 

After t = 115 ks (until the end of the observations), 
the count rates in the hard bands were gradually in- 
creasing, while the soft band (0.6-2 keV) did not show 
such a growth. This likely implies that the soft emis- 
sion has been heavily absorbed by dense matter located 
between the compact object and the observer. Large val- 
ues of the spectral hardness between the two softer bands 
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Fig. 4. — Light curves and hardness ratios from XIS around the 
flare at t = 110 ks. The time bin size is 28 s. Other conditions 
of the light curve extraction are the same as Figure \3\ At the two 
long data blanks from t = 106.8 ks and t = 111.8 ks, the satellite 
passed in the SAA, with the instruments turned off. Absence of 
data points that is seen between t = 107.4 ks and t = 111.8 ks is 
due to data removal by saturation of the satellite telemetry or by 
a condition of an elevation angle smaller than 5° above the Earth 
limb. 
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Fig. 5. — Light curves and hardness ratios of the XIS around 
the low state at t = 40 ks. Vela X-1 dropped into the low state 
from t = 38 ks until t = 41 ks. The time bin size is 28 s. Other 
conditions of the light curve extraction are the same as Figure |3] 



during that period also suggest strong absorption. Fast 
changes of this hardness ratio seen during that epoch 
can be caused by fluctuations of the absorbing matter 
density. 

3.3. Spectral Analysis 

In order to evaluate the spectral evolution of Vela X-1 
quantitatively, we have performed time-resolved spectral 
analysis. We extracted spectra and generated the de- 
tector response functions in the way described in ij3Aj 
Except for the soft excess due to the photoionized stellar 
wind, the X-ray emission below 10 keV is well described 
by an absorbed power law with iron fluorescence. Since 
in the preparatory analysis ( §3.ip wc have found the ab- 
sorption edge of iron at i^odgo ^7.1 keV, we model the 
soft X-ray emission by 

^ =exp(-7VH<Tabs(^))/cdgc(S)(AS-^+^Fc(£^)). (3) 

The additional factor fcd^c{E) to the preliminary model 
described by Equation ^ represents the effect of addi- 
tional absorption by iron, written as 



/cdgc(-B) — 



1 

CXp[-Tcdgc(-E/£;odgo) 



{E < Sedge) 
{E > Ecdgc), 



(4) 

where Todgo denotes the optical depth for the iron absorp- 
tion. 

We performed spectral fitting of the XIS spectra for 
intervals of 2268 s, which correspond to eight spin periods 
of the neutron star. Each spectrum has been extracted 
from the region of 45" < r < 120" to avoid the photon 
pileup effect. In the fitting, we have assumed iVjj, Tedge, 
the power-law index F, power-law normalization A, and 
the Fe line normalization to be free parameters. The 
energies of the iron line center and the iron K-edge were 
fixed to 6.4 keV and 7.25 keV, respectively, as have been 
determined by preliminary fitting. The energy range for 
fitting has been restricted to 2.5-10 keV to avoid the 
impact of the soft excess in the low energy band, as we 
mentioned in IJSHJ 

The results of the spectral fitting are tabulated in Ta- 
ble [21 and their time evolutions are plotted in Figure [Sj 
Reduced x^ values for each fit are also shown to check 
goodness of the fit. We excluded a few spectra with ex- 
posure time shorter than 568 s from the fitting analysis 
to avoid large uncertainties due to low statistics and po- 
tential bias by spin-phase dependent effects. This quan- 
titative analysis confirms the behaviors of the X-ray radi- 
ation from Vela X-1 seen in the light curves, as discussed 
previously. 

The X-ray ffux displayed strong variability while the 
power law index, F, remained relatively stable, with 
its value close to F ^ 1. Typically during the ob- 
servation, the flux ranged between 1 x 10~^ and 3 x 
10~^erg s^^cm^^, corresponding to ^ 100 mCrab in this 
band. As seen in the light curve (Fig. [3|), the object 
dropped into a very low-luminosity state with a flux of 
1.8 X 10~^°erg s~^cm~^ at < = 40 ks, showing a signifl- 
cant spectral softening to F = 2.2. In contrast, a sudden 
hardening to F = 0.5 occurred at t = 110 ks. Relatively 
hard states with F ~ 0.7 also appeared at t = ks and 
around t — 140 ks, when the flux was high. 
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The equivalent hydrogen column density -/Vjj estimated 
by X-ray attenuation is a good probe of the circumstel- 
lar matter in the Vela X-1 system. During the obser- 
vation, which corresponds to an orbital phase range of 
0.175-0.310, the obtained A'h value was showing strong 
variability between 1 x 10^^ and 9 x 10^^ cm^^. Since the 
interstellar absorption to Vela X-1 can be estimated as 
iVn ~ 6 X 10^^ cm~^ (assuming the distance of D = 1.9 
kpc and a typical hydrogen density of 1 cm~^), the time 
variability of A^h supports inhomogeneity of the absorb- 
ing matter inside the binary system. The strong stellar 
wind of the companion star, the accretion flow to the 
neutron star, and the accretion wake behind the neutron 
star can give significant contributions to the intrinsic ab- 
sorption in the system. 

The information about the circumstcllar matter inside 
the system is also reflected in the iron Ka line at 6.4 
keV, since the line is formed due to the reprocessing of 
the hard X-ray emission by matter. Roughly, the line 
intensity is proportional to the illuminating hard X-rays 
from the neutron star, and the line equivalent width is 
suitable for extraction of the information on the illumi- 
nated environment around the compact object. Though 
the equivalent widths have relatively large uncertainties, 
they seem to be correlated with A'h- The equivalent 
width became large, approaching a value of ^ 150 eV, 
around the epoch of t = 110 ks, when the spectrum shows 
the hardening, implying an increase of the matter density 
near the compact object. 

Since the absorption hardens the spectral shape, the 
estimated values of the photon index T and the column 
density iVu may be strongly coupled. Therefore, the 
compensating behaviors of A'h and T might be a result 
of wrong fitting due to the coupling. To study this we 
have selected the significant hardening at i = 110 ks (ID 
48) and the significant softening at i = 40 ks (ID 17), 
the two most prominent features detected for the photon 
index. For these two cases, we examined the x^ struc- 
tures in the two-parameter (F-A^h) space, and show the 
error contours for these fits in Figure [T] Although the 
two parameters have a correlation, the 99% confidence 
intervals of them are significantly smaller than the tem- 
poral changes of them at the significant hardening and 
softening. Thus, we conclude that the significant spectral 
variations are real, i.e., are not a result of the parameter 
coupling. 

We also checked correlations between spectral param- 
eters obtained over the whole observation. Figures [8] and 
[9] show the obtained relations between the X-ray flux and 
two parameters derived though the spectral fitting (iVu 
and the iron Ka line equivalent width). Both the X-ray 
absorption measured by A^h and the equivalent width of 
the iron Ka line do not show any significant correlations 
with the X-ray flux. By contrast, there is a positive cor- 
relation between Nn and the equivalent width of iron Ka 
line, as shown in Figure [101 

4. XIS/HXD WIDE-BAND ANALYSIS 

In this section we present the broad-band analysis of 
the data obtained with both XIS and HXD onboard 
Suzaku from the Vela X-1 pulsar. The key difference of 
this analysis to the one conducted in the previous section 
is that the spectral fitting was performed for a broader 
range of energies. Namely, the data obtained in the soft 



(2.5-10 keV) and hard X-ray (15-90 keV) energy bands 
were used simultaneously for the spectral fitting. The 
soft X-ray data were the same as in the previous section 
and the hard X-ray data were obtained with the HXD. 

4.1. Light Curves 

Figure [TT] shows light curves of the HXD-PIN data 
above --12 keV together with the XIS data (5-10 
keV). The curves were extracted from the two detectors 
with a time bin width of 284 s, which corresponds to 
the spin period (the same as Figure [3|). The HXD-PIN 
light curves are background-subtracted and dead-time- 
corrected though the non-Xray background of the HXD- 
PIN detector was ^0.5 counts s~^, which is not signifi- 
cant when compared to such a bright source as Vela X-1. 

While the temporal behavior in the two bands is very 
similar on the whole, the hardness ratio increases with 
the luminosity. The low states (low-luminosity states 
with spectral softening) found in the soft band are also 
clearly seen in the wide-band analysis (e.g. see a strong 
dip in the hardness ratio around t ^ AQ ks). The flare at 
f = 110 ks also showed a significant hardening. 

4.2. Spectral Model 

Since there is no well-established physical model of the 
radiation from accreting X-ray pulsars, we introduce a 
phenomenological spectral model to quantify the Suzaku 
wide-band data. For the detailed time-resolved analysis, 
we constructed a spectral model template by using the 
time-averaged spectrum of Vela X-1 obtained with the 
XIS and HXD. The spectra for the simultaneous wide- 
band modeling were extracted from data of XIS 0, XIS 3, 
HXD-PIN, and HXD-GSO with common time intervals 
for these detectors. 

There are three well-known phenomenological spectral 
models to describe a power law with quasi-exponential 
cutoff at hig h energy, i.e., a spectrum t ypical for X-ray 
pulsars (e.g. iKrevkenbohm et al.l Il999f). A power law 
with a high-energy cut is one of them (jWhite et aI1ll983D , 
written as 



dN 
'dE 



AE-^ 



exp 



{E < E,) 
E-EA ,_ _, (5) 



E 



'/ 



{E > E,) 



where E^ and Ef are referred to as "cutoff energy" and 
"folding energy", respectively. This model converges to 
the "normal" cutoff power law (power law with an expo- 
nential cut) if Ec -^ 0. Although this model has been 
often used for spectral modeling of the emission of ac- 
creting neutron stars, its derivative is discontinuous at 
E = Ec, and this could introduce artificial features in 
the continuum. 

The second one is a power l aw with the so-called Fermi- 
Dirac cutoff (|Tanakal Il98l IKrevkenbohm et al.l [19991 . 
which renders a spectrum with a smooth turnover: 



dN 
'dE 



AE- 



1 



exp 



E- Er 



E 



f 



(6) 



+ 1 



This function is smooth at any point, and nicely fits a 
spectrum with a cutoff. Due to historical reasons, the 
parameters of this model have the same names as those 
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Fig. 6. — Best-fit parameters of the XIS data to power-law models as functions of time. The values are tabulated in Table |2] The error 
bars represent la uncertainties. 
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Fig. 7. — Error contours in the F-A'^h spaces for the spectra ID 
48 (left) and ID 17 (right), which showed significant hardening and 
softening, respectively. The three contours are for 68%, 90%, and 
99% confidence levels. 



of the high-energy cut modeL However, the aetual mean- 
ing of these parameters is very different, therefore they 
cannot be compared directly. 

The third is the negative and posi tive exponen- 
tial cutoff power law (NPEX) model (jMiharal [l995l: 
iMakishima et al.lll999l ). 



S^'-^.-- 



p+r2 



Fi > 0, Ta > 0. 



-I ) , (7) 

(8) 

This is a combination of negative and positive power laws 
with an exponential cutoff by a common folding energy. 
The combination of the two power laws with largely dif- 
ferent indices provide good fits for a wide variety of spec- 
tral shapes. If the positive power-law index is r2 = 2 as 



f^*-.-- 



A2E-')cxp[~^ 



(9) 
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ID 



TABLE 2 

Results of spectral fit of the XIS data to power-law models 



Nh 



[10^ 



Photon 

r 



indc 



Power- 



law 



Fc norm. 
[10^^ ph cm^ 



Fo EW 

[cV] 



Flux (2-10 koV) X /dof 

[10^^ org cm^^s"^] 






1134 


4.83 ±0.18 


0.15 ±0.02 


0.73 ±0.03 


0.231 ±0.014 


4.86 ±0.42 


81 ± 7 


3.742 ±0.017 


1.12 


1 


3402 


3.10 ± 0.16 


0.07 ±0.02 


1.13 ±0.03 


0.305 ±0.017 


1.46 ±0.33 


39 ±9 


2.590 ±0.010 


0.99 


2 


5671 


2.54 ±0.11 


0.03 ±0.01 


1.09 ±0.02 


0.367 ±0.014 


2.41 ±0.29 


50 ±6 


3.484 ±0.010 


0.99 


4 


10207 


2.91 ±0.23 


0.06 ±0.03 


1.29 ±0.04 


0.121 ±0.010 


0.60 ±0.16 


54 ± 15 


0.797 ±0.005 


1.05 


5 


12475 


3.33 ±0.26 


0.11 ±0.04 


1.30 ±0.05 


0.168 ±0.015 


0.72 ±0.24 


48 ± 16 


1.051 ±0.008 


0.97 


6 


14744 


2.43 ±0.25 


0.05 ±0.03 


1.14 ±0.05 


0.230 ±0.020 


1.52 ±0.41 


55 ± 15 


2.013 ±0.015 


1.03 


7 


17012 


2.64 ±0.14 


0.07 ±0.02 


1.20 ±0.03 


0.253 ±0.013 


1.28 ±0.23 


47 ±9 


1.961 ±0.008 


0.93 


8 


19280 


2.43 ±0.26 


0.03 ±0.04 


1.23 ±0.05 


0.308 ±0.029 


1.80 ±0.48 


57 ± 16 


2.321 ±0.019 


1.00 


9 


21548 


1.88 ±0.14 


0.04 ±0.02 


1.14 ±0.03 


0.304 ±0.016 


2.19 ±0.31 


59 ±8 


2.758 ±0.011 


0.96 


10 


23817 


1.82 ±0.15 


0.09 ±0.02 


1.06 ±0.03 


0.272 ±0.014 


2.08 ±0.32 


55 ±9 


2.782 ±0.012 


0.96 


11 


26085 


2.63 ±0.27 


0.08 ±0.04 


1.17 ±0.05 


0.228 ±0.022 


1.68 ±0.41 


64 ± 16 


1.879 ±0.017 


0.99 


12 


28353 


3.13 ±0.20 


0.10 ±0.03 


1.30 ±0.04 


0.154 ±0.011 


0.68 ±0.18 


50 ±14 


0.965 ±0.006 


1.00 


13 


30621 


3.74 ±0.36 


0.12 ±0.05 


1.35 ±0.07 


0.115 ±0.015 


0.41 ±0.21 


44 ±23 


0.638 ±0.008 


1.05 


14 


32890 


3.18 ±0.19 


0.07 ±0.03 


1.25 ±0.04 


0.251 ±0.017 


1.17 ±0.28 


48 ± 12 


1.731 ±0.009 


0.96 


15 


35158 


2.69 ±0.17 


0.08 ±0.02 


1.25 ±0.03 


0.248 ±0.015 


1.37 ±0.25 


57 ±10 


1.744 ±0.008 


1.00 


16 


37426 


3.77 ±0.33 


0.04 ± 0.05 


1.55 ± 0.07 


0.119 ±0.014 


0.42 ±0.15 


63 ±22 


0.478 ± 0.005 


0.96 


17 


39694 


5.70 ±0.45 


0.06 ±0.07 


2.22 ±0.09 


0.157 ±0.026 


0.09 ±0.08 


37 ±33 


0.183 ±0.003 


0.96 


18 


41963 


3.62 ±0.33 


0.13 ±0.05 


1.37 ±0.07 


0.124 ±0.015 


0.49 ±0.21 


51 ±22 


0.677 ±0.007 


0.89 


19 


44231 


4.04 ±0.22 


0.07 ±0.03 


1.35 ±0.04 


0.157 ±0.012 


0.74 ±0.18 


57 ± 14 


0.879 ±0.006 


0.98 


20 


46499 


2.87 ±0.20 


0.10 ±0.03 


1.26 ±0.04 


0.161 ±0.012 


1.06 ±0.20 


68 ± 14 


1.110 ±0.006 


0.91 


21 


48768 


3.05 ±0.22 


0.09 ±0.03 


1.20 ±0.04 


0.164 ±0.013 


1.35 ±0.24 


76 ± 13 


1.252 ±0.008 


0.96 


22 


51036 


2.81 ±0.16 


0.07 ±0.02 


1.18 ±0.03 


0.199 ±0.011 


1.20 ±0.21 


54 ±9 


1.594 ±0.007 


1.04 


23 


53304 


1.77 ±0.16 


0.02 ±0.02 


1.12 ±0.03 


0.248 ±0.014 


1.56 ±0.31 


50 ± 10 


2.359 ±0.011 


1.08 


24 


55572 


1.39 ±0.13 


0.03 ±0.02 


1.12 ±0.03 


0.260 ±0.012 


1.89 ±0.26 


58 ±8 


2.498 ±0.010 


0.96 


25 


57841 


2.37 ±0.17 


0.08 ±0.02 


1.32 ±0.03 


0.215 ±0.013 


1.21 ±0.20 


66 ± 11 


1.373 ±0.006 


0.98 


26 


60109 


1.66 ±0.18 


0.03 ±0.03 


1.26 ±0.04 


0.260 ±0.017 


1.78 ±0.29 


71 ± 12 


1.949 ±0.010 


1.02 


27 


62377 


1.98 ±0.15 


0.04 ±0.02 


1.24 ±0.03 


0.219 ±0.012 


1.66 ±0.21 


76 ±9 


1.649 ±0.007 


0.98 


28 


64645 


2.27 ±0.22 


0.04 ±0.03 


1.29 ±0.04 


0.163 ±0.013 


0.84 ±0.21 


56 ±14 


1.109 ±0.007 


0.98 


29 


66914 


3.08 ±0.24 


0.02 ±0.03 


1.48 ± 0.05 


0.169 ±0.015 


0.84 ±0.18 


77 ±16 


0.808 ±0.006 


1.00 


30 


69182 


2.16 ±0.27 


0.08 ± 0.04 


1.29 ±0.05 


0.114 ± 0.011 


0.74 ±0.19 


71 ± 19 


0.783 ±0.007 


1.05 


31 


71450 


3.55 ±0.32 


0.00 ±0.09 


1.49 ±0.06 


0.179 ±0.020 


0.62 ±0.23 


55 ±22 


0.809 ±0.010 


0.96 


32 


73718 


3.62 ±0.23 


0.01 ±0.04 


1.54 ±0.05 


0.181 ±0.015 


0.55 ±0.16 


53 ± 15 


0.750 ±0.006 


0.97 


34 


78255 


5.52 ±0.26 


0.05 ± 0.03 


1.43 ±0.05 


0.197 ±0.017 


1.04 ±0.21 


75 ± 16 


0.893 ±0.006 


1.01 


35 


80523 


3.50 ±0.17 


0.04 ± 0.02 


1.30 ±0.03 


0.252 ±0.015 


1.57 ±0.23 


70 ± 10 


1.579 ±0.007 


1.01 


37 


85060 


3.39 ±0.18 


0.04 ± 0.02 


1.32 ±0.03 


0.207 ± 0.013 


0.99 ±0.20 


55 ± 11 


1.274 ±0.006 


1.07 


38 


87328 


4.92 ±0.26 


0.08 ±0.03 


1.29 ±0.05 


0.212 ±0.019 


1.23 ±0.28 


64 ± 15 


1.249 ±0.009 


1.05 


39 


89596 


4.74 ±0.27 


0.06 ±0.03 


1.21 ±0.05 


0.154 ±0.014 


1.22 ±0.25 


74 ± 16 


1.074 ±0.008 


0.96 


40 


91864 


4.66 ±0.20 


0.07 ±0.03 


1.26 ±0.04 


0.182 ±0.012 


1.22 ±0.20 


70 ±12 


1.150 ±0.006 


0.99 


42 


96401 


4.32 ±0.21 


0.10 ±0.03 


1.19 ±0.04 


0.148 ±0.011 


1.12 ±0.20 


69 ± 12 


1.072 ±0.006 


0.98 


43 


98669 


5.74 ±0.26 


0.12 ±0.03 


1.32 ±0.05 


0.172 ±0.015 


1.21 ±0.21 


81 ±14 


0.935 ±0.007 


0.98 


44 


100937 


4.71 ±0.32 


0.09 ± 0.04 


1.14 ±0.06 


0.205 ±0.022 


1.53 ±0.42 


62 ± 19 


1.602 ±0.015 


1.02 


45 


103206 


4.25 ±0.15 


0.09 ±0.02 


1.15 ±0.03 


0.249 ±0.013 


2.26 ±0.24 


76 ±9 


1.962 ±0.008 


0.98 


46 


105474 


4.01 ±0.23 


0.09 ±0.03 


1.14 ±0.04 


0.222 ±0.017 


2.35 ±0.33 


87 ± 12 


1.805 ±0.011 


0.92 


47 


107742 


8.23 ±0.27 


0.16 ±0.03 


0.97 ±0.04 


0.149 ±0.012 


2.54 ±0.28 


104 ± 12 


1.372 ±0.009 


1.05 


48 


110010 


4.13 ±0.22 


0.23 ±0.02 


0.46 ±0.04 


0.084 ± 0.006 


5.36 ±0.31 


149 ±9 


2.272 ±0.012 


1.07 


49 


112279 


6.59 ±0.31 


0.10 ±0.03 


1.01 ±0.05 


0.211 ±0.021 


3.18 ±0.43 


99 ±14 


1.939 ±0.016 


0.98 


50 


114547 


6.53 ±0.18 


0.08 ±0.02 


1.14 ±0.03 


0.221 ±0.013 


2.47 ±0.23 


92 ±9 


1.628 ±0.008 


1.03 


51 


116815 


4.91 ±0.25 


0.13 ±0.03 


1.16 ±0.05 


0.161 ±0.013 


1.58 ±0.25 


84 ±13 


1.195 ±0.008 


1.05 


52 


119083 


4.64 ±0.22 


0.15 ±0.03 


1.05 ±0.04 


0.143 ±0.011 


1.74 ±0.24 


85 ± 11 


1.313 ±0.007 


0.96 


53 


121352 


4.91 ±0.18 


0.14 ±0.02 


1.16 ±0.03 


0.209 ±0.013 


2.01 ±0.23 


83 ± 10 


1.542 ±0.007 


0.99 


54 


123620 


4.61 ±0.21 


0.15 ±0.03 


1.10 ±0.04 


0.189 ±0.014 


2.25 ±0.27 


92 ± 11 


1.575 ±0.009 


0.95 


55 


125888 


7.03 ±0.19 


0.11 ±0.02 


1.05 ±0.03 


0.196 ±0.012 


3.04 ±0.23 


110 ±9 


1.633 ±0.007 


1.04 


56 


128156 


7.47 ±0.22 


0.11 ±0.02 


1.08 ±0.04 


0.213 ±0.015 


2.76 ±0.27 


96 ± 10 


1.676 ±0.009 


0.96 


57 


130425 


4.93 ±0.16 


0.17 ±0.02 


0.93 ±0.03 


0.201 ±0.011 


3.51 ±0.26 


99 ± 7 


2.227 ±0.009 


1.03 


58 


132693 


6.76 ±0.16 


0.17 ±0.02 


0.82 ±0.03 


0.193 ±0.010 


4.52 ±0.26 


107 ±6 


2.471 ±0.009 


0.92 


59 


134961 


8.30 ±0.22 


0.19 ±0.02 


0.79 ±0.04 


0.179 ±0.012 


4.80 ±0.32 


116 ±8 


2.306 ±0.011 


1.02 


60 


137229 


7.52 ±0.17 


0.19 ±0.02 


0.76 ±0.03 


0.166 ±0.009 


4.38 ±0.26 


107 ±7 


2.325 ±0.009 


0.90 


61 


139498 


9.20 ±0.21 


0.25 ±0.02 


0.63 ±0.03 


0.135 ±0.008 


5.44 ±0.27 


131 ±7 


2.220 ±0.009 


1.09 


62 


141766 


5.50 ±0.16 


0.14 ±0.02 


0.85 ±0.03 


0.256 ± 0.013 


5.82 ±0.32 


110 ±6 


3.249 ±0.012 


1.00 


63 


144034 


2.78 ±0.13 


0.14 ±0.02 


0.87 ±0.03 


0.224 ±0.010 


4.52 ±0.28 


101 ±6 


3.069 ± 0.011 


0.93 



All errors represent la uncertainties, a) normalization of the power law in units of photons kcV~ cm~ s~ at 1 keV. b) normalization of Fc Kc 
at 6.4 keV. c) equivalent width (EW) of Fe Ka at 6.4 kcV to underlying continuum. The data set referred in the text via the ID numbers arc 
marked with gray color. 
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Fig. 8. — Relation between A^h and the X-ray flux. 
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Fig. 10. — Relation between the equivalent width of iron Ka line 
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Fig. 9. — Relation between the equivalent width of iron Ka line 
and the X-ray flux. 
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Fig. 11.— Light curves of XIS (5-10 keV) and HXD-PIN above 
~ 12 keV. A hardness ratio between the different detectors are also 
presented. The time bin size is 284 s, the same as the spin period 
of the neutron star. The HXD-PIN light curve are background- 
subtracted and dead-time-corrected. 



the function is known to be a good approxima- 
tion of the unsaturated Co mptoniza ti on sp ectrum 
(|Sunvaev fc TitarchukI [l980l) . IMiharal ([l99l found 
that X-ray pulsar spectra obtained with the Ginga 
sateUite are well des cribed by this model (sec also 
iMakishima et al.lll999| ). The NPEX model approaches 
to a normal cutoff power law when A2 — >■ 0. 

Examining these three models, we select the NPEX 
function (Eq. [9|) to characterize the spectrum and to ex- 
tract physical information. The high-energy cut model 
(Eq. [5]) is excluded because of the discontinuity of 
the derivative. Although the Fermi-Dirac cutoff model 
(Eq. [6]) well describes the spectral shape, results of pre- 
liminary fitting to the obtained data have shown that 
there is strong coupling between the power-law normal- 
ization A and the cutoff energy E^ near E'c ~ 0, which 
generates large uncertainties of these parameters. The 
NPEX function is smooth and easy to analyze because it 
simply consists of two power laws with a common ex- 



ponential cutoff. Furthermore, the model parameters 
allow a tentative physical i nterpretation in the context 
of thermal Comptoniza tion (jSunvaev fc TitarchukI II 980l : 
IMakishima et al.lll999l ). 

Many X-ray pulsars display absorption-like features 
in the hard X-ray continuum called cyclotron reso- 
nance scattering features (CRSFs). They are a re- 
sult of scattering via cyclotron resonance of electrons 
whose energy is quantized to the Landau levels in 
the strong magnetic field {B - lO^^ G). The CRSFs 
have complex structures which are composed of ab- 
sorption and emission via cyclotro n radiation from 
electrons in hot magnetized plasma (lArava fc Harding 
19991: lArava-Gochez fc Hardi^l 120001 : iSchonherr et al 



20071 ). It is, however, difficult to include these com- 



plex models of the CRSF in the spectral fitting with 
the detector response. We therefore adopt a widely- 
used phenomenolo gical model for cyclotron absorp- 
tions (e.g. .Mihara et all 119901: IMakishima et al.l 119991: 



10 
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iKrevkenbohm et alj 120081) . The effective optical depth 
due to the cyclotron absorptions can be written as 



TcycliE) =Di 



(E - S,yel)2 + Wi^ 
{W2E/2E,y,xY 



(10) 



£). 



[E - 2E, 



cyclj 



where Ecyd is the cyclotron energy, Wi, W2, I?i, D2 de- 
notes widths and depths for the fundamental (subscript 
1) and the first (subscript 2) harmonic resonances. The 
effect of the absorption in the spectral fitting is then in- 
cluded by the factor exp(— Tcyci). 

Combining the continuum model with the effects of 
absorption and emission from the circumstellar medium 
such as the stellar wind which is described by Equa- 
tion ([3|), we introduce a spectral model, based on the 
NPEX function, for the wide-band fitting as the follow- 



ing: 



dN 
'dE 



: exp(-iVHCrabs(-B))/odgo(^) 



{AiE 



-r 



A2E')cxp[ -^] +Ffc{E) 



Ef 



X exp(-Tcyci(£^)). 



(11) 



Since the model described by Equation (jlip is very com- 
plex and the x^ minimization might converge to an un- 
desired local minimum, we preliminarily fitted the global 
structure of the spectrum, and then set the initial values 
of the parameters of subsequent fits to the parameters 
obtained in the previous fit. 

The fitting result of the time-averaged spectrum to 
the NPEX model is shown in Table [H First we tried 
to fit the spectrum to the model without the CRSFs, 
but the data-to-model ratio showed a dip at £■ ~ 50 
keV, which is probably a CRSF at 50-55 keV t hat has 
been reported in many observations fe.g.lOr landi n et al.l 
199a IMakishima et al.l 119991: IKrevkenbohm et al.l 119991 
20021 ). In addition, some obser vations reported a.n ex- 



istence of a CRSF at 25 keV (jMakishima et al.l 119991: 
IKrevkenbohm et al.l 120021) while the observation with 
Beppo-SAX could not find a strong evidence of t he fea- 
ture around E ^ 25 keV (jOrlandini et al.l 119981 ). Our 
preliminary fit also showed a shallow dip at 25 keV. Con- 
sidering these facts, we added two CRSFs with the cy- 
clotron energy Ecyd — 25 keV (equal to the fundamen- 
tal line) to the preliminary model. The final results are 
shown in Figure [12] The model also contains the Fe K/3 
line at i? = 7.1 keV, but the significance of this line is 
not high. 

There are still residuals between the model and the 
data, probably, because the adopted spectral model, 
which is a combination of the NPEX function and the 
simple absorption-like functions describing the CRSFs. 
Indeed it is obvious that this simple model is not enough 
to represent the complicated integration of the neutron 
star emission over the long time and over the different 
emission regions. 

We used a spectral model with two CRSFs, the fun- 
damental line at ^^25 keV and the first harmonic line at 
-50 kcV, for fitting the data. While the CRSF at £; = 50 



keV is clearly seen in the spectrum, we do not find any 
strong evidence of the CRSF at 25 keV for the time- 
averaged spectrum. Thus, the shallow dip at 25 keV, 
which was attributed to a CRSF at the preliminary anal- 
ysis, is likely to related a discrepancy between the simple 
continuum model and the cor nplicated neutron sta r emis - 
sion. However, we note that IKrevkenbohm et al.l (|2002l ) 
showed the existence of the fundamental line at 25 keV 
based on spin-phase-resolved spectroscopy, since this fea- 
ture is much weaker than the CRSF at 50 keV and its 
visibility seems to depends on the pulsar phase. Thus, 
to confirm or discard its existence it would be impor- 
tant to perform the phase-resolved analysis, but such a 
study requires an additional dedicated analysis, and will 
be discussed elsewhere. 

4.3. Time-Resolved Spectra 

Since the hard X-ray radiation of Vela X-1 displays a 
strong time variability, it is important to analyze wide- 
band X-ray spectra time-resolved on a short time scale. 
We extracted wide-band spectra from the event data of 
XIS 0, XIS 3, XIS 1, and HXD-PIN for time bins of 2268 
seconds (the same approach as for the XIS analysis de- 
scribed in i i3.3p . In order to characterize the time evolu- 
tion of the neutron star radiation, we performed spectral 
fitting of the obtained spectra to the NPEX model de- 
scribed by Equation (|lip . Due to the lower statistics 
of the time-resolved spectra as compared to the time- 
averaged one, the parameters related to the CRSF were 
fixed to the values obtained from the fit of the time- 
averaged spectrum. 

The results of the fitting are tabulated in Table |3] and 
plotted in Figure [131 We excluded several spectra with 
exposure times shorter than 568 s from the fitting anal- 
ysis to avoid large uncertainties due to low statistics and 
potential bias by the pulsar spin-phase. Each obtained 
spectrum is well represented by the NPEX model. The 
energy fluxes in 2-10 keV and iVn estimated here are 
consistent with the values obtained only from the soft 
X-ray data (Table [2]). Though the absolute values of the 
photon index of the NPEX model cannot be compared 
directly with the index of a single power law, the tem- 
poral behavior is very similar to the one obtained in the 
soft-X-ray analysis. 

Together with the photon index F, the normalization 
of the positive power law A2 and the cutoff energy Ef 
characterize the hard X-ray radiation of the neutron star. 
The photon index F determines the slope of the flat 
power-law spectrum at lower energies. Figure [TJ] shows 
the relation between the photon index and the flux, and 
the shown data imply that the spectrum becomes hard 
when the flux increases. 

The two parameters of the NPEX model, A2 and Ef 
determine the spectral shape of the quasi-exponential 
cutoff at high energies, and therefore should reflect the 
physical conditions of the accreted plasma that radiates 
hard X-rays. The relation between A2 and the flux is 
shown in Figure [151 If the spectral shape were not chang- 
ing with the total luminosity, A2 would be exactly pro- 
portional to the flux level. But, interestingly, the positive 
power-law component is small and does not increase with 
flux if the flux is smaller than -^ 1.5 x 10~^ erg cm~^s~^. 
At higher fluxes, A2 increases with the flux, like the "nor- 
mal" negative power-law component. This behavior of 



TABLE 3 

Results of spectral fit of the time-averaged wide-band spectrum to the NPEX model 
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Parameter 

A^H [10^^ cm-' 

£,dgc [keV] 



""V5K 



Negative power- law photon index i' 

Negative power-law normalization Ai [photons s~^cm~^ at 1 keV] 

Positive power-law normalization A2 [10~^ photons s~^cm~^ at 1 keV] 

Cutoff energy Ef [keV] 



2.495 ± 0.043 
7.255 ±0.017 
0.095 ± 0.004 



0.385 ±0.016 

0.1076 ±0.0019 

0.152 ±0.013 

7.10 ±0.21 



Fe Ka energy [keV] 
Fe Ka width a [keV] 
Fe Ka normalization [10~^ photons s~^c 
Fe K/9 energy [keV] 
Fe K/3 width a [keV] 
Fe K/J normalization [10^'^ photons s~^c 



6.410 ±0.008 
0.00 ± 0.0001 
1.847 ±0.025 
7.08 ±0.16 
0.00 ±0.04 
0.109 ±0.025 



Fundamental cyclotron line depth Di 

Fundamental cyclotron line energy -Bcycl [keV] 

Fundamental cyclotron line width Wi [keV] 

First harmonic cyclotron line depth D2 

First harmonic cyclotron line energy 2_E(.ycl [keV] 

First harmonic cyclotron line width W2 



0.036 ± 0.008 
24.67 ±0.25 

9.0 ±1.4 

0.528 ±0.057 

49.34 ± 0.50 (fixed to 2£cycl) 

9.0 ±1.4 (fixed to Wi) 
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Fig. 12. — Time-averaged wide-band X-ray spectrum of Vela X-1 and the NPEX model fitted to the data. The data of XIS (black), 
XIS 3 (red), HXD-PIN (green), and HXD-GSO (blue) are shown. The top panel shows unfolded uF,^ spectrum with the detector response 
(top) and the bottom shows the data-to-modol ratio. The dashed linos show the two power-law components of the NPEX function. The 
model includes the CRSFs. 



the positive power law is apparent in the eorrclation plot 
between A2 and Ai , as shown in Figure 1161 

The cutoff energy Ef determines the shape of the high- 
energy turnover. The estimated values of Ef range be- 
tween 6 keV and 10 keV, and are concentrated around 
6-7 keV. The relation between the cutoff energy and 
the flux is shown in Figure [TT] It seems to be weakly 
antieorrelated with the flux if the flux is below ~ 2 x 
10^^ erg cm^^s^^. For higher fluxes, the cutoff energy 



seems to be constant about Ef ^ Q keV. 

5. DISCUSSIONS 

X-ray absorption and reprocessing provide us with cru- 
cial information on circumstellar medium in an HMXB. 
Since the surroundings of the neutron star in the Vela X- 
1 system is an ideal site to investigate fueling processes to 
the accreting system as well as feedback processes of the 
accretion power, the circumstellar medium has been one 



12 



Odaka et al. 



9. 


CM 


3 


o 


1 




', 




? 




Q) 





6 
4 
2 

1.5 
1 

0.5 




10 




6 - 



0.2 - 



0.1 



§ g 0.001 
-.^ o 

iSs: 0.0005 






20 



40 



60 



80 



100 



120 



140 



time [ks] 



Fig. 13. — Best-fit parameters of the XIS/HXD data to NPEX models as functions of time. The values are tabulated in Table|4] 
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Fig. 15. — Relation between the normalization of the positive 
power law A2 and the X-ray flux. 
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TABLE 4 
Results of spectral fit of the XIS/HXD data to NPEX models 



Time 



Nh 



[10^ 



Photon index ' 

r 



Cutoff energy 
Ef [kcV] 



Negative- PL 
norm. Ai 



Positive-PL 
norm ^2'^ 



Flux (2-10 
[10~ erg cm 



kcV) 

-2„-li 



x"/dof 






1134 


5.00 ±0.31 


0.15 ±0.07 


6.23 ±0.10 


0.17 ±0.02 


1.02 ±0.09 


3.713 ±0.027 


1.15 


1 


3402 


2.10 ±0.20 


0.39 ±0.04 


6.91 ±0.26 


0.18 ±0.01 


0.18 ±0.04 


2.597 ±0.014 


0.99 


2 


5671 


1.37 ±0.12 


0.32 ±0.02 


7.13 ±0.18 


0.20 ±0.01 


0.21 ±0.03 


3.513 ±0.011 


1.03 


4 


10207 


2.24 ±0.25 


0.65 ±0.05 


7.54 ±0.38 


0.08 ±0.01 


0.04 ±0.01 


0.781 ±0.007 


1.07 


5 


12475 


2.14 ±0.27 


0.55 ± 0.05 


7.60 ±0.56 


0.09 ±0.01 


0.05 ± 0.02 


1.074 ± 0.011 


1.00 


7 


17012 


1.99 ±0.18 


0.56 ±0.04 


6.28 ±0.10 


0.18 ±0.01 


0.30 ±0.02 


1.994 ±0.008 


1.01 


9 


21548 


1.41 ±0.19 


0.56 ± 0.04 


6.46 ±0.11 


0.24 ±0.02 


0.42 ±0.03 


2.846 ±0.013 


1.01 


10 


23817 


1.22 ±0.18 


0.44 ±0.04 


6.46 ±0.12 


0.20 ±0.01 


0.38 ±0.03 


2.819 ±0.013 


1.07 


11 


26085 


2.01 ±0.35 


0.52 ± 0.07 


6.88 ±0.28 


0.15 ±0.02 


0.19 ±0.04 


1.821 ±0.020 


0.98 


12 


28353 


2.16 ±0.21 


0.61 ±0.04 


7.83 ±0.37 


0.09 ±0.01 


0.05 ±0.01 


0.982 ±0.006 


1.11 


13 


30621 


3.07 ±0.43 


0.79 ±0.08 


8.07 ±0.99 


0.08 ±0.01 


0.02 ±0.01 


0.577 ±0.013 


1.03 


14 


32890 


2.43 ±0.23 


0.58 ± 0.04 


7.18 ±0.26 


0.16 ±0.01 


0.12 ±0.02 


1.736 ±0.011 


0.99 


15 


35158 


1.65 ±0.18 


0.53 ±0.03 


7.54 ±0.29 


0.14 ±0.01 


0.09 ±0.02 


1.775 ±0.009 


1.03 


16 


37426 


2.70 ±0.36 


0.81 ±0.07 


7.58 ± 0.68 


0.07 ±0.01 


0.02 ±0.01 


0.490 ±0.008 


1.08 


17 


39694 


5.48 ±0.46 


1.81 ±0.09 


11.22 ± 1.85 


0.13 ±0.02 


0.00 ±0.00 


0.186 ±0.005 


0.87 


18 


41963 


3.34 ±0.50 


0.78 ±0.10 


6.75 ±0.48 


0.09 ±0.02 


0.05 ± 0.02 


0.608 ±0.013 


0.98 


19 


44231 


3.36 ±0.28 


0.68 ±0.05 


7.66 ±0.41 


0.10 ±0.01 


0.04 ±0.01 


0.870 ±0.008 


1.00 


20 


46499 


1.87 ±0.22 


0.56 ±0.04 


7.79 ±0.39 


0.09 ±0.01 


0.05 ±0.01 


1.129 ±0.007 


0.96 


21 


48768 


2.50 ±0.32 


0.55 ±0.07 


6.69 ±0.23 


0.12 ±0.02 


0.14 ±0.02 


1.325 ± 0.014 


1.00 


22 


51036 


2.07 ±0.25 


0.52 ± 0.05 


6.21 ±0.14 


0.15 ±0.01 


0.25 ±0.03 


1.703 ±0.011 


0.98 


23 


53304 


1.56 ±0.27 


0.57 ±0.06 


6.16 ± 0.14 


0.20 ± 0.02 


0.37 ±0.04 


2.259 ±0.017 


1.07 


24 


55572 


0.87 ±0.19 


0.53 ± 0.04 


6.34 ±0.10 


0.21 ±0.02 


0.42 ±0.03 


2.713 ±0.012 


1.00 


25 


57841 


1.74 ±0.23 


0.69 ± 0.05 


6.92 ±0.18 


0.14 ±0.01 


0.12 ±0.01 


1.349 ±0.008 


0.97 


26 


60109 


1.28 ±0.29 


0.62 ±0.06 


6.44 ±0.20 


0.18 ±0.02 


0.21 ±0.03 


1.906 ±0.017 


1.17 


27 


62377 


1.01 ±0.22 


0.53 ± 0.05 


6.78 ±0.20 


0.15 ±0.01 


0.17 ±0.02 


1.837 ±0.011 


0.99 


28 


64645 


1.24 ±0.26 


0.57 ± 0.05 


8.10 ±0.68 


0.09 ±0.01 


0.04 ±0.02 


1.165 ±0.013 


1.08 


29 


66914 


2.06 ±0.26 


0.75 ±0.05 


7.38 ±0.40 


0.10 ±0.01 


0.04 ±0.01 


0.825 ± 0.007 


0.97 


31 


71450 


2.78 ± 0.50 


0.77 ±0.10 


6.46 ±0.61 


0.11 ±0.02 


0.05 ±0.02 


0.760 ±0.019 


1.11 


32 


73718 


2.55 ±0.19 


0.80 ± 0.04 


7.66 ±0.65 


0.10 ±0.01 


0.02 ±0.01 


0.762 ±0.007 


1.02 


34 


78255 


4.72 ±0.22 


0.84 ±0.11 


9.98 ±2.42 


0.12 ±0.01 


0.01 ±0.01 


0.893 ±0.032 


1.05 


35 


80523 


2.55 ±0.19 


0.60 ± 0.03 


7.76 ±0.39 


0.14 ±0.01 


0.06 ±0.02 


1.547 ±0.009 


1.02 


37 


85060 


2.46 ±0.15 


0.67 ±0.03 


8.91 ±0.79 


0.13 ±0.01 


0.03 ±0.01 


1.296 ±0.009 


0.96 


38 


87328 


4.53 ±0.28 


0.63 ± 0.05 


8.06 ±0.84 


0.13 ±0.01 


0.04 ±0.02 


1.202 ±0.017 


1.09 


39 


89596 


3.97 ±0.28 


0.61 ± 0.05 


8.72 ± 1.25 


0.10 ±0.01 


0.03 ±0.02 


1.064 ± 0.018 


0.96 


40 


91864 


3.41 ±0.20 


0.50 ± 0.03 


7.75 ±0.46 


0.10 ±0.01 


0.05 ±0.01 


1.171 ±0.008 


0.98 


42 


96401 


3.30 ±0.19 


0.51 ± 0.05 


8.64 ± 1.11 


0.08 ±0.01 


0.02 ±0.02 


1.073 ±0.015 


0.93 


43 


98669 


4.44 ±0.20 


0.61 ±0.08 


9.42 ± 1.79 


0.09 ±0.01 


0.01 ±0.01 


0.954 ±0.020 


0.98 


45 


103206 


3.13 ±0.17 


0.39 ±0.03 


7.03 ±0.21 


0.14 ±0.01 


0.14 ±0.02 


1.998 ±0.008 


1.01 


46 


105474 


2.89 ±0.29 


0.39 ±0.06 


6.91 ±0.27 


0.13 ±0.01 


0.17 ±0.03 


1.884 ±0.016 


1.01 


47 


107742 


7.42 ±0.32 


0.27 ±0.05 


7.28 ±0.37 


0.09 ±0.01 


0.11 ±0.03 


1.384 ±0.013 


1.10 


48 


110010 


2.87 ±0.29 


-0.31 ±0.06 


7.14 ±0.39 


0.04 ±0.00 


0.20 ±0.09 


2.137 ±0.020 


1.10 


49 


112279 


5.94 ±0.42 


0.39 ±0.09 


6.42 ±0.18 


0.15 ±0.02 


0.38 ± 0.05 


2.005 ± 0.026 


1.11 


50 


114547 


5.97 ±0.24 


0.53 ± 0.05 


6.38 ±0.10 


0.17 ±0.01 


0.30 ±0.02 


1.655 ±0.008 


1.00 


51 


116815 


3.95 ± 0.24 


0.46 ±0.04 


7.86 ± 0.60 


0.10 ±0.01 


0.05 ±0.02 


1.215 ±0.012 


1.04 


52 


119083 


3.59 ±0.23 


0.31 ±0.04 


7.84 ±0.64 


0.08 ±0.01 


0.05 ±0.02 


1.315 ±0.011 


1.15 


53 


121352 


3.75 ±0.20 


0.39 ±0.04 


6.98 ±0.23 


0.12 ±0.01 


0.13 ±0.02 


1.575 ±0.008 


0.92 


54 


123620 


3.45 ±0.25 


0.33 ±0.05 


6.77 ±0.26 


0.10 ±0.01 


0.15 ±0.03 


1.595 ±0.011 


1.01 


55 


125888 


6.06 ±0.24 


0.34 ± 0.05 


6.57 ±0.14 


0.12 ±0.01 


0.24 ± 0.03 


1.660 ±0.009 


1.03 


56 


128156 


6.17 ±0.33 


0.42 ±0.07 


6.38 ±0.15 


0.15 ±0.02 


0.32 ±0.04 


1.794 ±0.015 


1.06 


57 


130425 


4.22 ±0.26 


0.32 ±0.06 


6.43 ±0.12 


0.14 ±0.01 


0.41 ±0.04 


2.259 ±0.014 


1.12 


58 


132693 


5.90 ±0.23 


0.18 ±0.05 


6.21 ±0.08 


0.13 ±0.01 


0.61 ±0.04 


2.498 ±0.010 


1.03 


59 


134961 


7.67 ±0.38 


0.12 ±0.07 


7.03 ±0.23 


0.11 ±0.01 


0.31 ±0.06 


2.316 ±0.023 


1.09 


60 


137229 


5.82 ±0.27 


0.15 ±0.06 


6.43 ±0.12 


0.12 ±0.01 


0.51 ± 0.05 


2.482 ±0.015 


1.09 


61 


139498 


7.32 ±0.33 


0.08 ±0.07 


6.45 ±0.14 


0.11 ±0.01 


0.52 ± 0.07 


2.442 ±0.018 


1.04 


62 


141766 


5.67 ±0.25 


0.34 ± 0.05 


6.34 ±0.09 


0.22 ±0.02 


0.77 ±0.05 


3.268 ±0.017 


1.10 


63 


144034 


2.60 ±0.21 


0.30 ± 0.05 


6.15 ±0.08 


0.17 ±0.02 


0.68 ±0.05 


3.007 ±0.014 


1.02 



All errors represent la uncertainties, a) photon index of the negative power law (PL), b) normalization of the negative power law in units of 
photons cm~ s~ at 1 keV. c) normalization of the positive power law in units of 10~' photons keV~ cm~ s~ at 1 kcV. The data set referred in 
the text via the ID numbers are marked with gray color. 
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Fig. 16. — Relation between the normalizations of the two power 
laws, A2 and Ai. 
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Fig. 17. — Relation between the normalization of the cutoff en- 
ergy Ef and the X-ray flux. 



of the main s ubjec t for observations o f Vela X-1 (e.g. 
Nagase et all 119861 : iSako et all 119991 : iWatanabe et al.l 
20061 ). Based on the time variabihty of the absorption 
in Vela X-1, the relation between the structure of the 
stellar wind and the variability of the accreting object 
is discussed in this section. In addition, we briefly dis- 
cuss the radiation from the central engine based on the 
wide-band X-ray spectral data. 



5.1. Time Variability of Wind-fed Accretion System 

Typically, X-ray flux levels of wind-fed HMXBs show 
high variability, which is also seen in the light curves of 
Vela X-1 obtained with Suzaku (see Figs. El [TT|) . Simi- 
lar flaring behavior in hard X-rays was already observed 
with INTEGRAL up to an inten sity of ~ 5 Crab at 20-40 
keV (jKrevkenbohm et al.l 120081) . Since hard X-rays also 
display similar variability to that of the soft band, the 
variability should be caused by changes of the intrinsic 



radiation from the neutron star, not due to an increase 
of the attenuation by the surrounding material around 
the neutron star. 

Inhomogencity of the stellar wind is a possible origin 
of the time variability (Sec Section 1). The luminosity 
Lx by the wind-fed accretion varies as 



^X OC Pwind ^rol 



(12) 



where Pwind and v^ei are the density and the rela- 
tive velocity of the wind near the compact object 
(|Bondi fc HovldfToIl . 

It is probable that the wind of an OB star has highly 
clumpy structure in which density fluctuates with a fac- 
tor of 3-10 originated fr om the line-driven ins t abilit y 
of stellar UV photons (Pcssart fc Owockil [20031 |2005[ ). 
The porosity of the wind is also inferred from sym- 
metrical shapes of X-ray emission lines from isolated 
0-type stars obtained by high-resolution spe ctroscopy 
(|Oskinova et al.l [20061 : lOwocki fc CoIIct] ^QM) ■ Stellar 
winds in HMXBs are more complex due to the presence 
of the compact object and the strong X-ray radiation. 
ASCA data of Vela X-1 in its eclipse phase suggested 
that the stellar wind of Vela X-1 consists of a number 
of dense clumps ernbedde d in a tenuous, highly ion ized 
plasma (jSako et al.lll999[ ). IWatanabe et alJ (|2006( ) did 
not find any strong evidence for the clumpy structure of 
the wind in Chandra high- resolution data of Vela X-1, 
but confirmed that the wind velocity is suppressed near 
the neutron star due to the photoionization by the pow- 
erful X-rays, suggesting a strong impact of the X-rays on 
the dynamical structure of the wind. 

We examine here the Suzaku data based on the wind- 
clumping hypothesis. We would expect some absorption 
features in soft X-ray bands associated with increase of 
the mass accretion rate if a dense clump dropped into 
the neutron star. However, there is no significant trend 
of increase of X-ray absorption accompanied by luminous 
phases (Fig. [6]). In Figure [181 we also show relation be- 
tween A'h and the calculated X-ray luminosity Lx both 
of which are based on the XIS/HXD simultaneous anal- 
ysis. There is no apparent correlation between the two 
quantities, and moreover TVh can be a very small value of 
(1 — 2) X 10^^ cm~^ even at high-luminosity states with 
Ljf - 5 X 10^6 erg s'^. 

If the short-time variability of the fiux and circumstel- 
lar absorption is a reflection of the wind clumpy struc- 
ture, it is possible to estimate the physical properties of 
the individual clumps. Indeed, the light curves obtained 
with Suzaku (see Fig. [3]) suggest the flare duration tf 
of 1-10 ks. Assuming the relative velocity of the wind 
I'rci — 400 km s~^ at the neutron star, we obtain the 
clump radius of R^ ~ Vrcit//2 = (2 - 20) x 10^° cm. 
Since this radius is smaller than the accretion radius of 
Vela X-1 Race ~ 3 x 10^^ cm, the neutron star can accrete 
the whole clump. If tf = 1 ks and the flare luminosity 
Lx = 10"^^ erg s~^ are assumed, the total energy Lxtf 
relates to the clump mass Mc, as 



Lxt 



xif 



GMMc 



Mr_^ 



LxtfRtf 



= 4x10 



19 , 



(13) 



where Af, and i?, are the mass and the radius of the 
compact star, respectively. Thus, we obtain the mean 
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Fig. 18. — Relation between the equivalent hydrogen column den- 
sity A^H a-iid the X-ray luminosity of Vela X- 1 . The two quantities 
are calcula ted based on the XIS/HXD simultaneous analysis de- 
scribed in ^4.31 The luminosity is integrated energy over 0.1-100 
keV, but the data fitting is based on 2.5-60 keV. 



number density of the clump as 
A'L 3 



TOp AttRc^ 



7 X 10" cm"3 



(14) 



and the radial column density Nc of the clump 1.4 x 
10^^ cm~^. Consequently, the mean number density and 
the column density can be written as 



ric =7 X 10 



11 



Lx 



1 ks 



J^Q37 Qj-g g 1 
-2 



Wrcl 



400 km s' 



Nc =ncRc 



1 ks 



1.4 X 10 

-1 



22 



Lx 



I'rcl 



1037 erg s-i 

-2 



(15) 



(16) 



400 km s- 



Since the number density of the wind around the neu- 
tron star which is estimated by the standard wind model 
(jCastor. Abbott, fc KleiiJ [l975l ) is n = 3 x 10^ cm-^, 
the clump density estimated by equation (|15p for a short 
flare oi tf '^ I ks is higher by two orders of magnitude. 
However, hy drodynamical simulations by iBlondin et al.l 
(|1990l I1991D show existence of density fluctuation near 
the compact object with a factor of '-^100, and therefore 
such a dense clump may cause a short flare. The column 
density of the clump is enough to detect by observation 
when the flare is bright and short. The luminous flare 
at i = 110 ks (see ij3.2p is a candidate for such a flare. 
We roughly estimate the column density of the clump 
that would generate the flare. Adopting the luminosity 
Lx = 10^7 erg s~^ and the flare duration tf — 400 s esti- 
mated by visual inspection of the light curve (Fig. [4]), we 
obtain the number density ric = 4 x lO'^^ cm"^ and the 
column density Nc = 3.5 x 10^^ cm~^. The difference of 
observed column densities just before the flare (ID: 47) 
and during the flare (ID: 48) is AA'^h = 4 x 10^^ cm~^. 
The good agreement of the decrease of the column den- 
sity with the estimated value of the clump supports that 
the clump responsible for the absorption causes the fol- 



lowing flare. 

A larger wind clump can generate a longer flare. We 
observed several moderate flares which have luminosi- 
ties of Lx ~ 5 x 10^^ erg s~^ and durations of i/ ~ 5 
ks. Such a flare corresponds to the number density 
He = 1.4 X 10^° cm~3 and the column density Nc = 
1.4 X 10^^ cm~^. Since this density is five times as dense 
as the value of the smooth wind model, the density en- 
hancement is naturally explained by the line-driven in- 
stability of the stellar wind. The column density of the 
clump, however, is too small to detect by X-ray observa- 
tions because of contamination of other absorption effects 
due to complicated wind structure close to the neutron 
star, as discussed in the next subsection. 

In summary, the Suzaku observation of Vela X-1 shows 
consistency with the wind clumping responsible for the 
time variability of the X-ray radiation. If the short 
flares of Vela X-1 are triggered by such clumps, they 
are likely to be very dense clumps around the neutron 
star, not the clumps which possibly exist in wide re- 
gions of the stellar wind. However, other mechanisms 
such as the instability of the temporary accretion disk 
(iTaam fc Frvxell 1989) or th e so-called "flip-flop insta- 
bility" (jMatsuda et al.iri991| ) can possibly also explain 
the flaring behavi or of the wind-fed accretin g pulsar (See 
also discussion in iKrevkenbohm et al.|[2008l ). 

5.2. Enhanced Wind 

It has been well known that Vela X-1 shows s trong vari- 
ability of absorption with the orbital phase (jPan et al.l 
119941 : iNagase et al.l |198"6[ ). Chandra observations also 
showed very different column densities of TVh = 1.45 x 
10^2 cm-2 at phase 0.25 and jVh = 18.5 x 10^^ cm^^ 
at phase 0.5 (jWatanabe et al.l 120061 ). These absorp- 
tion features are attributed to enhanced winds or ac- 
cretion wakes around the neutron star, which are gener- 
ated by competing effects including gravitational, rota- 
tional, and radiation press ure forces, and X-ray heating 
(jBlondin et al.l 119901 I1991D . Suzaku allows us to trace 
the structure of such dense clouds close to the neutron 
star with finer time and spatial resolutions than those of 
the past observations. 

The Suzaku observation covered the orbital phase be- 
tween 0.175 and 0.31. This early orbital phase has been 
thought to be a period in which the X-ray absorption 
is small, as observed with the Chandra at phase 0.25. 
However, the light curves and the time-resolved spectral 
analysis of the Suzaku data shows relatively heavy ab- 
sorption in the second half of the observation. Similar 
absorption increase were seen by t he Tenma satellite a t 
the phase <j) = 0.25 (See Fig. 2 of INagase et~aIlll986D . 
The Suzaku results also show strong fluctuation of ab- 
sorption after t ~ 100 ks (i.e. at orbital phase 0.27). 

As a possible interpretation, we attribute this absorp- 
tion feature to a bow shock generated by the interaction 
of the moving compact object with the su personic stellar 
wind. The hydrodynamic simulations by IBlondin et al.l 
(|1991| ) show such a structure emerging around the com- 
pact object. According to the simulation, dense regions 
at the downstream of the shock will be observed as the 
heavy absorption at orbital phase (j) ~ 0.4. The related 
hump of absorption is in fact very similar to the absorp- 
tion increase observed by Suzaku (See Figures 2 and 8 
from .Blondin et al., ,19911) . Although the orbital phase 
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of the predicted and observed features are different, tliis 
interpretation still remains feasible, since the position of 
the bow shock greatly depends on the physical conditions 
of the binary system, including the binary separation, the 
wind properties, and so on. 

5.3. Spectral Features of the Central Engine 

As shown in Figure [TH we have found a correlation 
between the X-ray luminosity and the spectral hardness. 
Instead of the X-ray flux. Figure [!§] shows the same re- 
lation plotted with the luminosity in 2-10 keV of the 
source. The source spectrum becomes harder as the lu- 
minosity increases. On the assumption of the thermal 
Comptonization as the dominant process of the X-ray 
radiation, this correlation could have a natural interpre- 
tation: an enhanced accretion rate should lead both to an 
increase of the X-ray luminosity and to a more efficient 
Comptonization occurring in the accreted plasma. 

The spectral features of Vela X-1 — the power law with 
a quasi-exponential cutoff and the hardening with the 
luminosity — agree with the Comptonization hypothesis. 
In order to confirm it, it is necessary to obtain more quan- 
titative physical relations based on a more detailed emis- 
sion model. Such a model should incorporate a dynam- 
ical structure of the accretion flow and a magnetic field 
with physically r easonable values. In our subsequently 
pubhshed paper (|Odaka et al.ll2013[) . detailed modeling 
using Monte Carlo simulations and implications of the 
observation based on the modeling are presented. 




4 6 

Lx[10^'*ergs-i] 

Fig. 19. — Relation between the photon index F and the X-ray 
luminosity. Two significantly hard and soft points are out of the 
y-range. 



6. CONCLUSIONS 

We studied the brightest wind-fed accreting neutron 
star Vela X-1 with the Suzaku X-ray observatory. Due 
to the wide-band coverage and the low-noise performance 
of Suzaku, we have successfully obtained the time vari- 
ability of the wide-band X-ray spectrum over a short 
timescale of 2 ks for the observation duration of 145 ks. 



To exploit the high performance of Suzaku in the soft 
X-ray energy band, we have performed the data analysis 
in two steps: first we have analyzed the data obtained 
with XIS only; and after that we conducted a broad-band 
analysis, which accounted for both XIS and HXD data. 
The soft X-ray analysis has shown that: 

• The "low states", in which the source luminosity 
becomes an order of magnitude lower than the aver- 
aged value, still shows X-ray pulsations, indicating 
that the accretion flow reaches the magnetic poles 
of the neutron star. In this state, the source dis- 
plays a significant spectral softening with a photon 
index of ^2. 

• The circumstellar absorption seems not to correlate 
with the X-ray luminosity except for a very hard 
flare which occurred at i = 110 ks. This is con- 
sistent with the clumpy wind scenario to explain 
the time variability. However, other mechanisms 
such as the magnetospheric barrier effects cannot 
be excluded at this stage. 

• The equivalent hydrogen column density TVh in- 
creased in the second half of the Suzaku observa- 
tion, which corresponds to the orbital phases of 
0.25-0.3. This feature can be interpreted as a dense 
region where a bow shock is generated by interac- 
tions of the compact object with the stellar wind. 

The broad-band X-ray analysis allowed us to obtain 
the following conclusions: 

• The time-averaged wide-band spectrum is well rep- 
resented by the NPEX function, which is a com- 
bination of two different-slope power laws with a 
common exponential cutoff, with an obvious CRSF 
at 50 keV. However, we could not find any signif- 
icant evidence for a CRSF at 25 keV in this spin- 
phase-averaged spectrum. 

• The time- resolved wide-band spectra between 2.5 
keV and 50 keV are well described by the NPEX 
function. The time variability of the spectral prop- 
erties should reflect variability of the physical con- 
ditions of the accreted plasma. 

• Data obtained from Vela X-1 show correlation of 
the photon index F with the X-ray luminosity. On 
the assumption of the thermal Comptonization as 
the dominant process of the X-ray radiation, this 
correlation agrees with a natural physical relation 
that the optical thickness of the accreted plasma 
increases with the mass accretion rate onto the neu- 
tron star. 
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